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Abstract 

Increased evidence of porcine deltacoronavirus (PDCoV) causing diarrhoea in pigs 
has been reported in several countries worldwide. The virus has currently evolved 
into three separated groups including US, China and Southeast Asia (SEA) groups. In 
Vietnam, PDCoV was first reported in 2015. Based on phylogenetic analyses of spike, 
membrane and nucleocapsid genes, it is suggested that Vietnam PDCoV is chimeric 
virus. In the present study, we retrospectively investigated the presence of PDCoV in 
Vietnam and the full-length genomes of six PDCoV isolates identified in 2014-2016 
were further characterized. The results demonstrated that Vietnam PDCoV was first 
detected as early as 2014. All six Vietnam PDCoV are in the SEA group and further 
divided into two separated subgroups including SEA-1 and SEA-2. Vietnam PDCoV 
in SEA-2 was closely related to Thai and Lao PDCoV. Recombination analysis demon¬ 
strated that three isolates in SEA-1 were a chimeric virus of which P12_14_VN_0814, 
the first Vietnam isolate, and US PDCoV isolates were major and minor parents, re¬ 
spectively. The recombination was further evaluated by phylogenetic construction 
based on 3 recombinant fragments. The first and third fragments, closely related 
to P12_14_VN_0814, were associated with ORFla/lb and N genes, respectively. 
The second fragment, associated with S, E, and M genes, was closely related to US 
PDCoV isolates. High antigenic and hydrophobic variations were detected in SI pro¬ 
tein. Three-day-old pigs challenged with the chimeric virus displayed clinical diseases 
and villus atrophy. In conclusion, Vietnam PDCoV is genetically diverse influenced by 
an external introduction from neighbouring countries. The chimeric Vietnam PDCoV 
can induce a disease similar to Thai PDCoV. 
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1 | INTRODUCTION 

Porcine deltacoronavirus (PDCoV) is an enveloped, positive-sense, 
single-stranded RNA virus belonged to the genus Deltacoronavirus, 
family Coronaviridae, order Nidovirales (Woo, Huang, Lau, & Yuen, 


2010). PDCoV is the causative agent of a newly emerged enteric 
disease in pigs characterized by watery diarrhoea, dehydration and 
death with 40%-80% mortality following an outbreak (Jung, Hu, & 
Saif, 2016). PDCoV infects enterocytes of the small intestine lead¬ 
ing to villous atrophy and malabsorption (Hu et al., 2015; Jung et 
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al., 2015). Clinical signs caused by PDCoV resemble those with two 
viruses in the genus Alphacoronavirus, family Coronaviridae, includ¬ 
ing porcine epidemic diarrhoea virus (PEDV) and transmissible gas¬ 
troenteritis virus (TGEV) (Jung et al., 2015), but severity of clinical 
disease is lesser and piglet mortality is lower. PEDV and TGEV infect 
pigs of all ages and cause 100% morbidity, with 80%-100% mortality 
rates in piglets. The higher severity of PEDV and TGEV infection in 
pigs indicates that the two viruses are more suitable for infecting 
pigs. On the other hand, PDCoV may not be perfectly adapted to 
pigs, and the virus may be evolving towards better ability to infect 
a specific host. Recently, an experimental study in calves suggested 
that calves are susceptible to infection by PDCoV but not PEDV 
(Jung, Hu, & Saif, 2017). 

The full-length genome of PDCoV is approximately 25 kb in 
length and contains seven essential genes, including open reading 
frame la/lb (ORFla/lb), spike (S), envelop (E), membrane (M), and 
nucleocapsid (N), and specific accessory genes including non-struc- 
tural protein 6 (Nsp6) and Nsp7, flanked by a 5'- and 3'-untranslated 
region (UTR) (Woo et al., 2012). ORFla and ORFlb, occupying two- 
thirds of the genome, encode 2 overlapping viral replicase/transcrip- 
tase polyproteins, la and lb, which are cleaved into 15 Nsps, namely 
Nsp2-16 (Zhang & Yoo, 2016). The S gene encodes S protein which 
consists of two domains called SI and S2 domain. The SI domain 
plays an important role in binding to specific host cell receptors and 
contains a neutralizing epitope. The S2 domain functions in mem¬ 
brane fusion (Shang et al., 2018). The E and M genes are transmem¬ 
brane proteins associated with viral envelope formation and virus 
release (Woo et al., 2010). On the other hand, the functions of the 
N gene are associated with viral RNA replication and pathogenesis 
(Lee & Lee, 2015). 

PDCoV was first detected in Hong Kong (HKU15-44 and 
HKU15-155 isolates) in 2012 (Woo et al., 2012). Although it was first 
detected, there was no association with clinical disease. The first ev¬ 
idence of PDCoV causing a disease was first evident in Ohio, USA 
in 2014. Soon after the first emergence in Ohio, USA, the virus was 
then reported in 18 states (Jung et al., 2016). The rapid widespread 
of PDCoV was evident when the virus was subsequently reported 
in China, South Korea, Thailand, Laos, and Vietnam in 2015 (Le et 
al., 2018; Lee & Lee, 2014; Madapong et al., 2016; Marthaler, Jiang, 
Collins, & Rossow, 2014; Saeng-Chuto, Lorsirigool, et al., 2017; Song 
et al., 2015; Wang, Byrum, & Zhang, 2014a, 2014b). Although the 


first detection of PDCoV in those countries was in 2014-2015, the 
retrospective investigation of intestinal samples demonstrated that 
the presence of PDCoV in China, USA and Thailand was as early as 
2004, 2013 and 2013, respectively (Saeng-Chuto, Stott, et al., 2017; 
Sinha, Gauger, Zhang, Yoon, & Harmon, 2015). 

Presently, PDCoV has evolved into three separated groups, in¬ 
cluding US, China and Southeast Asia (SEA). PDCoV isolates from 
Thailand and Laos were clustered in a novel PDCoV group, SEA. The 
SEA group is genetically distinct from US and China PDCoV groups 
(Madapong et al., 2016; Saeng-Chuto, Lorsirigool, et al., 2017; Saeng- 
Chuto, Stott, et al., 2017). In Vietnam, PDCoV was first detected 
in 2015 (Saeng-Chuto, Lorsirigool, et al., 2017). Genetic analyses 
based on S, M, and N genes demonstrated that Vietnam PDCoV iso¬ 
lates were genetically different from the SEA group (Saeng-Chuto, 
Lorsirigool, et al., 2017). The phylogenetic analyses based on S and 
M genes demonstrated that Vietnam PDCoV is clustered in the SEA 
group. However, phylogenetic analysis based on the N gene demon¬ 
strated contradictory results, since Vietnam PDCoV was clustered 
in the US group. The findings implied that Vietnam PDCoV could 
potentially be a recombinant virus. Whether or not this genotype 
of PDCoV isolates predominantly existing in Vietnam remains not 
known. We therefore conducted a retrospective study investigating 
the presence of PDCoV in intestinal samples collected from pigs less 
than 5 days old experiencing severe diarrhoea, in Vietnam from 2011 
to 2016. The full-length genome sequences of the detected Vietnam 
PDCoV were further characterized on the basis of heterogenicity, 
and recombination analyses together with virus infectivity in cell cul¬ 
ture and pathogenicity in piglets were performed. 

2 | MATERIALS AND METHODS 
2.1 | Samples and viruses 

One hundred and eight intestinal samples collected from 49 swine¬ 
herds with clinical diarrhoea outbreaks in different regions in Vietnam 
(Table 1) were collected for viral pathogens causing diarrhoea during 
2011-2016. All intestinal samples were individually minced using an 
Axygen™ Tissue Grinder (Thermo Fisher Scientific) and resuspended 
in IX phosphate-buffered saline (PBS; 0.1 M, pH 7.2). The suspen¬ 
sions were centrifuged at 10,000 x g for 10 min at 4°C, and the 
supernatant filtered through a 0.45-p.m nylon membrane (Corning 


Year 

Regions 

No. farms 

Total samples 

+ve PDCoV 

+ve PEDV 

+ve TGEV 

2016 

Middle 

2 

6 

0 

5 

0 


South 

2 

4 

3 

4 

0 

2015 

South 

14 

29 

6 

15 

0 

2014 

South 

21 

47 

2 

43 

0 

2013 

South 

3 

9 

0 

7 

0 

2012 

South 

3 

6 

0 

6 

0 

2011 

South 

4 

7 

0 

7 

0 

Total 


49 

108 

11 

87 

0 


TABLE 1 Intestinal sample number 
that was collected from different regions 
in Vietnam in 2011-2016 
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Inc.). Virus was isolated, and plaque purification was performed as 
described in Data SI. 

2.2 | Virus detection and full-length 
genome sequencing 

Viral RNA was extracted from the supernatant using Nucleospin 
Viral RNA Extraction Kit (Macherey-Nagel Inc.) and converted to 
cDNA using M-MuLV Reverse Transcriptase (New England Biolabs 
Inc.). All samples were tested with the specific primer for S gene of 
PEDV and specific primer for N gene of TGEV as in previous studies 
(Lee & Lee, 2014; Park et al., 2007; Song et al., 2015; Temeeyasen et 
al., 2014; Wang, Byrum, & Zhang, 2014a; Woo et al., 2010). PDCoV 
was detected by PCR using specific primer for M and N genes as pre¬ 
viously described (Wang, Byrum, & Zhang, 2014b). The PCR prod¬ 
ucts were electrophoresed at 100 V for 30 min on a 1% agarose gel, 
and the gel was stained with RedSafe™ nucleic acid staining solution 
(iNtRON Biotechnology Inc.) and examined under a UV light. The 
PCR products were purified using Nucleospin Gel and PCR Clean-up 
kit (Macherey-Nagel Inc.), and sequencing was performed by First 
BASE Laboratories Inc. PDCoV-positive samples were subjected to 
whole-genome sequencing using 26 primers (Table SI). All specific 
PCR products were purified using Nucleospin Gel and PCR Clean-up 
kit (Macherey-Nagel Inc.), and sequencing was done by First BASE 
Laboratories Inc., using an ABI Prism 3730XL DNA sequencer. 

2.3 | Phylogenetic and genetic analyses 


methods implemented in RDP4 software (Martin, Murrell, Golden, 
Khoosal, & Muhire, 2015). Major and minor parents and breakpoints 
also were detected, and recombination fragments were generated 
with the RDP4 software (Martin et al., 2015). 

Phylogenetic trees based on recombinant fragments were sep¬ 
arately constructed using the ML method with 1,000 bootstrap 
replicates, and substitution models were designated by the best-fit 
substitution model function implemented in MEGA software version 
6.0 (Tamura et al., 2013). TN93 + G+l model was used in phyloge¬ 
netic tree based on all recombination fragments construction. 

2.5 | Sliding window, antigenicity and 
hydrophobicity analyses 

The sliding window, hydrophobicity and antigenicity were evaluated 
among the chimeric Vietnam PDCoV isolates in this study, and sus¬ 
pected major and minor parents of the chimeric Vietnam PDCoV, 
respectively (Figure 2). 

The nt sequence variation sites were predicted using dnaSP 
version 6 with 100 bp and a step size of 25 bp (Rozas et al., 2017). 
Antigenicity and hydrophobicity were also analysed using IEDB 
Analysis Resource (http://tools.immuneepitope.org/main) with the 
Kolaskar and Tongaonkar method and ExPASy-ProtScale (http://web. 
expasy.org/protscale) with the Kyte-Doolittle method, respectively. 

2.6 | In vitro infection of recombinant 
Vietnam PDCoV 


To determine the genetic relationship between the six Vietnam 
PDCoV isolates and other PDCoV isolates, the nucleotide (nt) and 
amino acid (aa) sequences were aligned using the CLUSTALW pro¬ 
gram (Thompson, Higgins, & Gibson, 1994), together with 84 other 
PDCoV isolates available in GenBank (Table S2). The nt and aa iden¬ 
tities were determined using the sequence identity matrix function 
implemented in BioEdit software (Hall, 1999). The phylogenetic tree 
based on the full-length genome, open reading frame (ORF) la/lb, 
SEMN and S, M and N genes was separately built using the maximum 
likelihood (ML) method with 1,000 bootstrap replicates, and substi¬ 
tution models were selected by the best-fit substitution model func¬ 
tion implemented in MEGA software version 6.0 (Tamura, Stecher, 
Peterson, Filipski, & Kumar, 2013). TN93 + G+l model was used in 
phylogenetic trees based on full-length genome and SEMN gene 
construction. TN93 + G model was used in phylogenetic trees based 
on ORFla/lb and S genes construction. K2 + G and K2 + G+l models 
were used in the construction of phylogenetic trees based on M and 
N genes, respectively. 

2.4 | Recombination analysis and phylogenetic 
reconstruction based on recombination fragments 

Recombination of six Vietnam PDCoV isolates was further analysed 
together with the 84 reference PDCoV isolates using automated 
RDP, GENECONV, BOOTSCAN, MaxChi, CHIMAERA and SISCAN 


Viral growth curves of the three chimeric Vietnam PDCoV isolates 
and NT1_1215 isolate were compared. Each isolate was performed 
in quadruplicate. The NT1_1215 isolate is a Thai PDCoV isolated 
from a swineherd in Thailand experiencing a diarrhoea outbreak 
with piglet mortality at approximately 40%. The full-length genome 
is deposited in GenBank with accession number KX361345. 

The PDCoV-positive samples were subjected to virus isola¬ 
tion using swine testicular (ST) cells (ATCC CRL-1746). The viral 
isolation was approved by Institutional Biosafety Committee of 
Chulalongkorn University (CU-IBC) (IBC 1831062). In brief, ST cells 
were seeded in T25 flasks (Corning, NY, USA) containing Advanced 
Minimum Essential Medium (adv. MEM) supplemented with 10% 
heat-inactivated foetal bovine serum (FBS), 1% HEPES, 1% antibi¬ 
otic-antimycotic, 1% L-glutamine (Gibco) and 0.0001% pancreatin 
from porcine pancreas (Sigma-Aldrich) for each viral inoculation. 
After the cells reached 80% confluency, the cells were inoculated 
with each virus at a titre of 10 3 TCID 50 /ml. The infected cells were 
incubated at 37°C and 5% C0 2 . The viruses were harvested at dif¬ 
ferent time points (0, 12, 24, 36, 48, 60 and 72 hr post-inoculation 
(hpi.)) and stored at -80°C. 

Viral titre at each time point was determined. ST cells were cul¬ 
tured in 96-well plates (Corning, NY, USA) containing supplemented 
adv. MEM. Cells were inoculated with each virus at each time point 
at 10-fold serial dilutions (10°-10 7 ). Cytopathic effect was deter¬ 
mined daily, and the growth curves were plotted. 
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2.7 | In vivo infection of recombinant 
Vietnam PDCoV 

Eighteen 3-day-old healthy piglets were randomly divided into three 
groups: (i) mock control group (n = 6), (ii) chimeric Vietnam PDCoV 
(P29_15_VN_1215) isolate-infected group (n = 6), and (iii) NT1_1215 
isolate-infected group (n = 6). Each piglet was inoculated with 5 ml 
of each virus at a titre of 10 3 TCID 50 /ml. Clinical signs were observed 
daily. Three piglets per group were euthanized at 3 and 5 days post¬ 
inoculation (dpi.), and five segments of the small intestine, including 
duodenum, proximal jejunum, middle jejunum, distal jejunum and 
ileum, were collected and fixed in 10% formalin. Haematoxylin and 
eosin (H&E) staining was performed, and villus/crypt ratio was cal¬ 
culated. Comparison between each group was performed with one¬ 
way analysis of variance, followed by Tukey's multiple comparison 
test implemented in GraphPad Prism 7 (GraphPad Software Inc.), 
p < .05 was considered significant. 


3 


RESULTS 


3.1 | PDCoV detection 


One hundred and eight intestinal samples collected from pigs with 
diarrhoea were assayed for the presence of three viral pathogens 
including PEDV, TGEV and PDCoV using PCR assays (Table 1). Of 
108 intestinal samples tested, none were positive for TGEV. Eleven 
(10.19%) and 87 (80.56%) samples were positive for PDCoV and 
PEDV, respectively. All PDCoV-positive samples were also positive 
for PEDV. PDCoV was detected in Vietnam as early as 2014, but the 
most detection was mainly in 2015-2016. Of 11 samples, only six 
samples including P12_14_VN_0814, P29_15_VN_1215, P30_15_ 
VN_1215, P1_16_VN_0116, P19_16_VN_0416, and P20_16_ 

VN_0416 isolates were able to characterize their full-length genome 
sequences and further genetic analyses. 


3.2 | Full-length genome characterization 

The full-length genome sequence of the six Vietnam PDCoV iso¬ 
lates, recognized as P12_14_VN_0814, P29_15_VN_1215, P30_15_ 
VN_1215, P1_16_VN_0116, P19_16_VN_0416 and P20_16_VN_0416, 
were deposited in the GenBank database under accession numbers 
MH700628, KX998969, MH118333, MH118331, MH118332 and 
MH700629, respectively. The six isolates were varied in size, varying 
from 25,405-25,409 nt. The major difference is the deletion in nu¬ 
cleotides in ORFla/lb region. The full-length genome sequence of 
P12_14_VN_0814 isolate was 25,406 nt in length excluding poly(A) 
tail. P29_15_VN_1215, P30_15_VN_1215 and P1_16_VN_0116 iso¬ 
lates were 3 nt longer in which their full-length genome sequences 
were 25,409 nt in length excluding poly(A) tail. The full-length genome 
sequence of P19_16_VN_0416 and P20_16_VN_0416 isolates were 
25,405 nt in length excluding poly(A) tail. All isolates were character¬ 
ized by the same gene order of the 5'-UTR, ORFla/lb, S, E, M, Nsp6, 
N, Nsp7 and 3'-UTR. The length of each region is presented in Table S3. 


Pairwise nt and aa sequence identity values of each gene be¬ 
tween the six Vietnam PDCoV isolates together with two previously 
reported Vietnam PDCoV isolates, HaNoi6/2015 and Binh21/2015 
isolates, and other reference PDCoV isolates are presented in 
Table S4. Based on the full-length genome analysis, the P12_14_ 
VN_0814, P29_15_VN_1215, P30_15_VN_1215 and Pl_16_ 

VN_0116 isolates were closely related to previously reported 
HaNoi6/2015 and Binh21/2015 isolates (98.7%-100% nt) than 
Thai and Lao PDCoV isolates (98.1%-98.4% nt), US PDCoV isolates 
(97.5%-98.1% nt), Japan PDCoV isolates (97.2%-98.1% nt), South 
Korean PDCoV isolates (97.2%-98.0% nt) and China PDCoV iso¬ 
lates (97.1%-98.3% nt). In contrast, P19_16_VN_0416 and P20_16_ 
VN_0416 isolates were more closely related to the Thai and Lao 
PDCoV isolates (99.1%-99.8% nt identity) than the HaNoi6/2015 
and Binh21/2015 isolates (98.3%-98.4% nt), US PDCoV isolate 
(97.2%-97.5% nt), China PDCoV isolates (96.9%-97.7% nt), Japan 
PDCoV isolates (96.9%-97.3% nt) and South Korean PDCoV iso¬ 
lates (96.9%-97.3% nt). 

The ORFla/lb gene of all six Vietnam PDCoV isolates owns aa 
deletion and insertion in their genomes compared to other PDCoV 
isolates. The deletions of two ( 401 [_K 402 ) and three ( 758 PVG 760 ) aa 
were observed at positions 401 to 402 and 758 to 760, respectively. 
Interestingly, these 5 aa deletions in the six Vietnam PDCoV isolates 
were similar to other PDCoV isolates from SEA countries. The inser- 

C A 

tion of one aa ( N) at position 51 of the S gene was found in all six 
Vietnam PDCoV isolates similar to US PDCoV isolates. Based on the 
E, M and N genes, all six Vietnam PDCoV isolates did not show any 
aa deletion or insertion compared to other PDCoV isolates. 

3.3 | Phylogenetic analysis 

Phylogenetic trees constructed based on the full-length genome 
and ORFla/lb demonstrated that the six Vietnam PDCoV isolates 
were grouped in the SEA group together with Thai, Lao, and two 
previously reported Vietnam PDCoV isolates based on full-length 
genome and ORFla/lb gene sequences. The SEA group was further 
evolved into two subgroups including subgroup SEA-1 and SEA-2 
(Figure la,b). P19_16_VN_0416 and P20_16_VN_0416 together 
with Thai and Lao PDCoV isolates were clustered in the subgroup 
SEA-1. In contrast, P12_14_VN_0814, P29_15_VN_1215, P30_15_ 
VN_1215 and P1_16_VN_0116 were clustered in subgroup SEA-2 
(Figure la,b). Interestingly, PDCoV isolates in the subgroup SEA-2 
are from Vietnam only. 

Phylogenetic trees constructed based on SEMN and S genes 
demonstrated that P19_16_VN_0416 and P20_16_VN_0416, and 
P12_14_VN_0814 were clustered together with the SEA group, sub¬ 
group SEA-1 and SEA-2, respectively (Figure lc,d). However, three 
of the six Vietnam PDCoV isolates, including P29_15_VN_1215, 
P30_15_VN_1215 and P1_16_VN_0116, were clustered together 
with the US group (Figure lc,d). 

The phylogenetic tree based on N gene demonstrated that 
P12_14_VN_0814, P29_15_VN_1215, P30_15_VN_1215 and Pl_16_ 
VN_0116 were clustered together with the China group, differing 
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FIGURE 1 Phylogenetic trees based on nucleotide sequences of full-length genome (a), ORFla/lb (b), SEMN (c), S (d), M (e) and N (f) 
genes. The phylogenetic trees were constructed using the maximum likelihood (ML) method with 1,000 bootstrap replicates implemented 
in MEGA software version 6.0. Red dots represent the PDCoV isolates from Vietnam in this study. All phylogenetic trees, except M gene, 
consisted of 3 groups namely US, China, and Southeast Asia (SEA) groups. SEA group was divided into 2 subgroups namely SEA-1 and SEA-2 
subgroups 
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95% Breakpoint confidence interval 
99% Breakpoint confidence interval 
Tract of sequence with a recombinant origin 

P12_14_VN_0814 (Major Parent) - USA/Minnesota 159/2014 (Minor Parent) 

P1214 VN 0814 (Major Parent) - P2915VN1215, P3015VN1215, and P1_16_VN_0116 (Recombinant) 

— USA/Minnesota 159/2014 (Minor Parent) - P29_15_VN_1215, P30_l 5_VN_1215, and P1_16_VN_0116 (Recombinant) 

FIGURE 2 Recombination plot between P12_VN_14_0814, USA/Minnesotal59/2014 and three Vietnam PDCoV isolates: P29_15_ 
VN_1215, P30_15_VN_1215 and P1_16_VN_0116 isolates as the major parent, minor parent and recombinant, respectively 


from P19_16_VN_0416 and P20_16_VN_0416, which were clus¬ 
tered in the SEA group (Figure If). 

3.4 | Recombination analysis and phylogenetic 
reconstruction based on recombination fragments 

Three recombination fragments were generated. The first and 
second fragments are at nt 1-19,453 and 23,511-25,426 (second 
fragment) in which are locations of ORFla/lb and N genes, re¬ 
spectively. The third fragment is at nt 19,454-23,512 in which is 
located at S, E and M gene region (Figure 2). Recombination analy¬ 
sis suggested that three Vietnam PDCoV isolates including P29_15_ 
VN_1215, P30_15_VN_1215 and P1_16_VN_0116 isolates were 
recombination of the P12_14_VN_0814 isolate (major parent) at 


the first and third fragments, and USA/Minnesotal59/2014 isolate 
(minor parent) at the second fragment (Figure 2). The phylogenetic 
trees based on the first and third fragments demonstrated that 
these three chimeric Vietnam PDCoV isolates were closely related 
to the P12_14_VN_0814, Hanoi/2015 and Binh21/2015 isolates 
(Figure 3a,c). However, the phylogenetic trees based on the second 
fragment suggested that that the three chimeric Vietnam PDCoV 
isolates were closely related to US PDCoV isolates (Figure 3b). 

3.5 | Sliding window, hydrophobicity and 
antigenicity analyses 

The full-length genome sequences of the five isolates includ¬ 
ing P29_15_VN_1215, P30_15_VN_1215, P1_16_VN_0116, 


FIGURE 3 Reconstruction of phylogenetic trees based on the nucleotide sequence of three recombinant fragments, including first 
(a), second (b) and third (c) fragments. The phylogenetic trees were constructed using the maximum likelihood (ML) method with 1,000 
bootstrap replicates implemented in MEGA software version 6.0. Red dots represent the PDCoV isolates from Vietnam in this study. Blue 
triangle represents P12_VN_14_0814 isolate (major parent). Green triangle represents USA/Minnesotal59/2014 isolate (minor parent) 



























SAENG-CHUTO etal 


Transboundary and Emerging Diseases 


-Wiley 



(b) 


Lj' 

U 


LC260044_YMG/JPN/2014 
KR265850_USA/MlChigan448/2014 
KR265849_USA/Michigan447/2014 
- KJ481931_USA/lllinois121/2014 

I-KR265865_U S A/lowa459/2014 

1— KR150443_U$A/Arkansas61/2015 

-LC260045 HKD/JPN/2016 


- KM820765_KNU 14-04 
LC260041JWT/JPN/2014 
-1 C260038_AKT/JPN/2014 

- LC260043 OKN/JPN/2014 
- KR265857_USA/lllinois273/2014 
' KR265856_USA/lllinois272/2014 

- KJ584357 HKU15/KY4813 

- KR265864_USA/Mmnesota292/2014 
-KY364365_KN U16-07 


KR265861_USA/Nebraska210/2014 
r-^| KR265860_USA/Nebraska209/2014 
KJ584359_HKU15/NE3579 
KJ584355_HKU15/IL2768 
KR265852_USA/lllinois449/2014 
-KY926512 KNU16-11 


LC260042_MYZ/JPN/2014 

— KJ620016_l 1KU15/M 16148 
- KJ769231 OhioCVMI/2014 
KR265863_USA/Ohio445/2014 

L. KR265862JJSA/Ohio444/2014 
KJ584356_HKU 15/SD3424 
KR265848_USA/Mrnnesota214/2014 
KR265847 USA/Minnesota442/2014 
KR265851_USA/lndiana453/2014 
LC260040 GNM-2/JPN/2014 
J 1 LC260039_GNM-1/JPN/2014 
- KR265853_USA/Minnesota/2013 

KR265855_USA/Minnesota455/2014 
KR265854_USA/Minnesota454/2014 
KJ584358_HKU15/PA3148 

KR265858_USA/NorthCarolina452/2014 
— KX022602_PDCoV/USA/lowa136/2015 
KX022604_PDCoV/USA/Nebraska 137/2015 
KX022605_PDCoV/USA/Ncbraska145/2015 
KX022603_PDCoV/USA/Minnesota140/2015 

- ▼ KR265859 USA/Minnesota159/2014 
KY354363_DH1 
KY354364_DH2 

“ KX998969_P29_15_VN_1215 
MH118333_P30_1S_VN_1215 
MH118331_P1_16_VN_0116 

KY363868_CHN-GD16 05 
-MF095123 CHN-HG-2017 


■ KR 

T-t 
-1 

f 

- Kl 

k - ' 

KR: 

■ KJ584 

'Ll 


, KT266822_CH/Sichuan/S27/2012 
I KT266822_CH/Sichuan/S27/2012 

- MF280390_CHN-GD-2016 

KY065120_CHN/Tianjin/2016 




-KR131G21 CHJXNI2/2015 


KT021234_CH/SXD 1 /2015 




KP757892_CHN-JS-2014 

- MG242062_CHN-HeB1 -2017 
MF431743_SD 
MF041982_SH JS/SL/2016 

-KU665558 CHN-LYG-2014 

MF948005_HB-BD 
KU981059_NH 

KY293677 CH/JXJGS01/2016 


<1 




KT336560_CHN-HN-2014 
-MF431742 GD 


c 


JQ065043_HKU 15/HKU15-155 
JQ065042_HKU 15/HKU15-44 
KP757891_CHN HB 2014 
-KY513724_CH/Hunan/2014 

-KY513725_CH/Jiangsu/2014 


KX443143 CH 01 


KP757890 CHN-AH-2004 


M F642322_CHN/GS/2016/1 
I-MF642325 CHN/QH/2017/1 


■ MF642323 CHN/GS/2016/2 
1 -MF642324_CHN/GS/2017/1 


I KX834351_HaNoi6/2015 
A MH700628_P12_14_VN 0814 
KX834352_Binh21/2015 
KY363867_CHN-GD16-03 

KX361343_ST1_0213 

KX361344 ST2 0313 


i—kx; 

L, 


d 


KX361345 _NT1 1215 
MH700627_NT1 -0416/2016 
- KU984334_TT-1115 
KU051641_Swine/Thailand/S5011/2015 
KU051649_Swme/Thailand/S5015U2015 
y KX118627 P116BTL 
MH118332_P19_16_VN_0416 
MH700629 P20 16 VN 0416 


■rar 


TTST 

































































































































































































10 


Wiley- 


Transboundary and Emerging Diseases 


(c) 


KJ584358_HKU 15/PA3148 
> KR265858_USA/Nof1hCarolina452/2014 
KR265853_USA/Minnesota/2013 
KY364365_KNU 16 07 
-KY926512 KNU16-11 


LC260042_M YZ/JPN/2014 
I LC260040_GNM-2/JPN/2014 
_ 1LC260039J3NM-1/JPN/2014 
KJ584355_HKU15/IL2768 
KR265852_USA/lllinois449/2014 
KJ584359_HKU15/NE3579 
KR265861 USA/Nebraska210/2014 
KR265860_USA/Nebraska209/2014 
- KJ584357_HKU 15/KY4813 

— KX022604_PDCoV/USA/Nebraska137/2015 
KX022603_PDCoV/USA/Minne$ota140/2015 
- KX022605_PC)CoV/U S A/Nebraska145/2015 
KX022602_PDCoV/USA/lowa136/2015 

-KJ769231_OhioCVM1/2014 

r-4 KR265862_USA/Ohio444/2014 
- KR265863_USA/Ohio445/2014 
KR265851_USA/lndiana453/2014 
r- KJ584356 HKU15/SD3424 
71 ,KR265848_USA/Minnesota214/2014 
' KR265847_USA/Minnesota442/2014 
KR265850_USA/Michigan448/2014 
KR265849_USA/Michigan447/2014 
▼ KR265859_USA/Mir>nesota159/2014 
KR265864_USA/Minnesota292/2014 

-LC260045 HKD/JPN/2016 


LC260038_AKT/JPN/2014 
- LC260041JWT/JPN/2014 
L KM820765_KNU14 04 
KR265865_USA/lowa459/2014 
- KJ620016_HKU15/MI6148 
LC260043_OKN/J PN/2014 
— LC260044 YMG/JPN/2014 

H KY354363 DH1 
KY354364_DH2 

KR265857_USA/lllinois273/2014 
U KR265856 USA/lllinois272/2014 

— KJ481931_USA/Illinoisl 21 /2014 

KR150443_USA/Arkansas61/2015 
KR265855_USA/Minnesota455/2014 
KR265854_USA/Mir>nesota454/2014 
-KY363868 CHN-GD16-05 


KP757892_CHN JS-2014 
• KT266822 CH/Sichuan/S27/2012 
I KT266822_CH/Sichuan/S27/2012 
MF431743_SD 

““ 15 |—MG242062_CHN-HeB1-2017 
i— MF041982 SHJS/SL/2016 

— KR131621 CHJXNI2/2015 


— KU665558_CHN LYG 2014 

-KT021234 CH/SXO1/2015 

-KY513724 CH/Hunan/2014 


HT, 


MF280390_CHN-GD-2016 
KY065120_CHN/Tianjin/2016 
MF095123 CHN-HG-2017 




JQ065042_HKU15/HKU15-44 

KP757891_CHN-HB-2014 

KX834352_Binh21/2015 
r=j A M H700628_P12_14_VN_0814 
KX834351_HaNoi6/2015 
-KY363867_CHN GD16 03 


I MH118333 P30 15 VN 1215 
-44, ■ MH118331_P1_16_VN_0116 
I KX998969_P29_15_VN_1215 
I KX361343_ST1_0213 
1KX361344 ST2 0313 

-KU984334_TT-1115 

| KX361345_NT1_1215 
MH700627 NT 1 -0416/2016 

KU051641_Swine/Thailand/S5011/2015 
KU051649_Swine/Thailand/S5015U2015- 

- w KX118627_P 1_16_BTL 

M H118332 P1916VN0416 
MH 700629 P20 16 VN 0416 


KY513725_CH/Jiangsu/2014 
KX443143_CH-01 


IKU9 

H 


KU981059_NH 

KY293677 CH/JXJGS01/2016 


■ KT336560_CHN-HN-2014 
1-MF431742 GO 


- JQ065043_HKU15/HKU15-155 
— MF948005_HB-BD 


KP757890 CHN-AH-2004 


- MF642322 CHN/GS/2016/1 

■ MF642323_CHN/GS/2016/2 

MF642325_CHN/QH/2017/1 
MF642324 CHN/GS/2017/1 


"ii 


TW 


SAENG-CHUTO etal. 


P12_14_VN_0814 (major parent) and USA/Minnesotal59/2014 
(minor parent) were used in the sliding window analysis. Five sites, 
designated Rl, R2, R3, R4 and R5, demonstrated high degree of 
genetic variation between all five isolates. Three sites including Rl 
(nt position 1,326-1,425), R2 (nt position 2,101-2,225) and R3 (nt 
position 5,601-5,775) were located within ORFla/lb regions. Two 
sites including R4 (nt position 20,101-20,200) and R5 (nt position 
25,001-25,100) were located within the S region, and between the 
N and 3'-UTR regions (Figure 4). 

The aa of the S gene encoding the structural protein associated 
with binding and neutralizing epitopes were subjected for further 
antigenicity and hydrophobicity analyses. Six segments at positions 
aa 7-15, 36-44, 116-124,134-139, 346-352 and 546-554 demon¬ 
strated a difference in hydrophobicity and antigenicity indices be¬ 
tween the five isolates (Figure 5). 

3.6 | In vitro infection with recombinant 
Vietnam PDCoV 

Viral growth curves of the chimeric Vietnam PDCoV (P29_15_ 
VN_1215) and Thai (NT1_1215) isolates in ST cells were compared. 
The two isolates had a similar pattern of viral growth curves, but 
virus titres were different (Figure 6). Cytopathic effect, character¬ 
ized by cell rounding and clumping, with both PDCoV isolates were 
first observed at 12 hpi with a virus titre of 10 TCID 50 /ml (mean 
log 10 TCID 50 /ml ± SD; 1 ± 0). At 24 hpi, the titre of P29_15_VN_1215 
remained at 10 1 TCID 50 /ml (mean log 10 TCID 50 /ml ± SD; 1 ± 0), while 
thetitre of NT1_1215 increasedto 10 2 TCID 50 /ml(mean log 10 TCID 50 / 
ml ± SD; 2 ± 0). The titres of both PDCoV isolates slowly rose from 
24 to 48 hpi. At 36 and 48 hpi, the titres of P29_15_VN_1215 were 
10 2 (mean log 10 TCID 50 /ml ± SD; 2 ± 0.25) and 10 4 TCID 50 /ml (mean 
log 10 TCID 50 /ml ± SD; 4 ± 0.25), respectively, which were lower than 
those of NT1_1215 at 10 3 (mean log 10 TCID 50 /ml ± SD; 3 ± 0.25) and 
10 5 TCID 50 /ml (mean log 10 TCID 50 /ml ± SD; 5 ± 0.25), respectively. 
The titres of both PDCoV isolates rapidly declined to 10 1 TCID 50 / 
ml at 60 hpi (mean log 10 TCID 50 /ml ± SD; 1 ± 0.25). However, at 72 
hpi., the titre of NT1_1215 increased to 10 2 TCID 50 /ml (mean log 10 
TCID 50 /ml ± SD; 2 ± 0), but the titre of P29_15_VN_1215 remained 
at 10 1 TCID 50 /ml (mean log 10 TCID 50 /ml ± SD; 1 ± 0). 

3.7 | In vivo infection with recombinant 
Vietnam PDCoV 

Pigs in the mock control group of which orally inoculated with culture 
media displayed no clinical disease throughout the study. In contrast, 
pigs orally inoculated with either Vietnam PDCoV isolates (P29_15_ 
VN_1215) or Thai PDCoV isolate (NT1_1215 isolate) displayed simi¬ 
lar level of clinical diseases associated with PDCoV including severe 
diarrhoea, vomiting, dehydration, weakness and lethargy. Clinical 
diseases were observed in both Vietnam and Thai PDCoV infected 
groups at 1 dpi. However, no piglets died before euthanasia. 

At 3 and 5 dpi, 3 piglets in each group were necropsied 
and five parts of small intestine including duodenum, proximal 


FIGURE 3 (Continued) 
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FIGURE 4 Differences in full-length genome sequence between five isolates, including P29_15_VN_1215, P30_15_VN_1215, Pl_16_ 
VN_0116, P12_VN_14_0814 (major parent) and USA/Minnesotal59/2014 (minor parent), were analysed by sliding window method. The red 
arrows show high variation regions 
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FIGURE 5 Hydrophobicity and antigenicity differences were found in S protein between the five isolates studied: P29_15_VN_1215, 
P30_15_VN_1215, P1_16_VN_0116, P12_VN_14_0814 (major parent) and USA/Minnesotal59/2014 (minor parent) 


jejunum, middle jejunum, distal jejunum and ileum were col¬ 
lected for further histopathological examination. H&E staining 
was performed. All H&E-stained small intestine segments of 


piglets infected with P29_15_VN_1215 and NT1_1215 showed 
villous atrophy compared to that of the mock control group (Gl) 
(Figure 7). 
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FIGURE 6 Viral growth curves of the Thai PDCoV (NT1_1215) (green line) and the chimeric Vietnam PDCoV (P29_15_VN_1215) (red line) 
isolates. Statistical significance between groups annotated by same lower-case letters, p-value < 0.05 


Villus/crypt (V/C) ratios of piglets in each group were determined 
(Figure 8). At 3 dpi, the V/C ratio in the three jejunum segments of 
piglets infected with both PDCoV isolates was significantly lower 
compared to that of the mock control group (p < .05) (Figure 8b- 
d). At 5 dpi, the V/C ratio in middle and distal jejunum segments of 
piglets infected with both PDCoV isolates was significantly lower 
than that in piglets in the mock control group (p < .05) (Figure 8c,d). 
Although no significant difference was observed, the V/C ratio in all 
parts part in both infected groups trend to reduced when compared 
between 3 and 5 dpi within the groups (Figure 8a,c). 


4 


DISCUSSION 


In Vietnam, the first detection of PDCoV was reported in 2015 in the 
southern region of Vietnam (Saeng-Chuto, Lorsirigool, et al., 2017). 
Genetic analyses based on S, M and N genes discovered contradic¬ 
tory results suggesting that Vietnam PDCoV could be a recombinant 
virus. However, the recombination of the virus has not been con¬ 
firmed and the full-length genome was not attempt in the previous 
study. We therefore conducted a retrospective study investigating 
the presence of PDCoV in Vietnam, genetic diversity and character¬ 
ize the full-length genome of detected PDCoV. The isolated chimeric 
virus was further evaluated for its pathogenicity. 

In the present study, PDCoV was first detected in southern re¬ 
gion of Vietnam as early as 2014. The full-length genome sequences 
of six Vietnam PDCoV isolates were characterized. Phylogenetic 
analysis based on the full-length genome demonstrated that 
Vietnam PDCoV was genetically diverse. The six Vietnam PDCoV 
isolates were grouped in the SEA cluster together with the Vietnam 
(Hanoi/2015 and Binh21/2015), Thai and Lao PDCoV isolates which 
all have an amino acid deletion at positions 401 to 402 and 758 to 
760 located in the ORFla gene. However, the six Vietnam PDCoV 
isolates were separated into two sub-clusters recognized as SEA-1 


and SEA-2. Two of six isolates were clustered in SEA-1 together with 
Thai and Lao PDCoV isolates while other four isolates were clus¬ 
tered in SEA-2 together with Hanoi/2015 and Binh21/2015 isolates. 

Due to a separated origin, the findings from the present study 
suggest there were at least 2 possible external introductions 
of PDCoV into Vietnam. As we mentioned earlier, the Vietnam 
PDCoV isolated in 2016 (P19_16_VN_0416 and P20_16_VN_0416) 
shares high degree of genetic similarities and was clustered to¬ 
gether with Thai and Lao PDCoV in subgroup SEA-1, the virus 
isolates might be introduced from these neighbouring countries. 
Secondly, the external introduction of US-related PDCoV was 
also observed. This introduction seemed to influence the evo¬ 
lution of Vietnam PDCoV as evidenced by chimeric PDCoV in 
subgroup SEA-2. The Vietnam PDCoV isolates detected in 2014 
(P12_VN_14_0814) were clustered together with other Vietnam 
PDCoV isolates, HaNoi6/2015 and Binh21/2015, in subgroup 
SEA-2. The HaNoi6/2014 isolate was detected in Ha Noi province 
in the northern region of Vietnam since October 2015 and consid¬ 
ered as recombinant virus between Thai and China isolates (Le et 
al., 2018). This study shows the P12_VN_14_0814 isolate shares 
100% similarities with the HaNoi6/2014 isolate which means the 
P12_VN_14_0814 isolate could be also recombinant virus be¬ 
tween Thai and China isolates. Interestingly, three of other four 
isolates in subgroup SEA-2 were different. The P29_15_VN_1215, 
P30_15_VN_1215 and P1_16_VN_0116 isolates were confirmed 
to be chimeric virus which were a recombination between 
Vietnam PDCoV (P12_VN_14_0814), the major parent, and US 
PDCoV (USA/Minnesotal59/2014), the minor parent. As previous 
reports, the recombination in RNA virus occurs when a host cell 
is infected with at least two viral genomes at the same time and 
affects host cell tropism, immune evasion, virulence and patho¬ 
genesis (Martin et al., 2011; Simon-Loriere & Holmes, 2011). The 
findings in this study suggested that there are at least four differ¬ 
ent genome sequences of PDCoV currently circulating in swine 
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FIGURE 7 Haematoxylin and eosin (H&E) staining of five segments of small intestine, including duodenum (a-c), proximal jejunum (d-f), 
middle jejunum (g-i), distal jejunum (j-l) and ileum (m-o). Villous atrophy was observed in all small intestine segments of piglets infected with 
the chimeric Vietnam PDCoV (P29_15_VN_1215) (b, e, h, k and n) or the Thai PDCoV (NT1_1215) (c, f, i, I and o) isolates, compared to mock 
control group (a, d, g, j and m). Statistical significance between groups annotated by asterisk, p-value < 0.05 
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FIGURE 8 Villus to crypt (V/C) ratio of five segments of small intestine, including duodenum (a), proximal jejunum (b), middle jejunum (c), 
distal jejunum (d) and ileum (e) 


farms in Vietnam. The virus may be trying to evolve, adapt to pigs, 
emerge and cause a disease. It is remained to be seen whether 
each of four different genome sequence groups would influence 
the genetic evolution of each other and become a dominant group, 
or continuously recombine and evolve in the population, resulting 
in the emergence of a novel viral isolate. 

To analyse the nt variation between chimeric Vietnam PDCoV 
isolates and their parents, the sliding window method was used. The 
results showed that high nt variation was detected in five sites of the 
full-length genomes. Most nt variations occurred in the ORFla/lb 
gene. One of the five sites was found in the S gene which encodes 
the S protein. The hydrophobicity and antigenicity variations of the S 
protein were then analysed. Hydrophobicity is associated with pro¬ 
tein-protein interaction (e.g. ligand of virus-host receptor binding) 
and can be changed by replacement with hydrophobic amino acids 
(Bosch, Zee, Haan, & Rottier, 2003; Tekewe, Connors, Middelberg, 
& Lua, 2016), and change in antigenicity can affect immune re¬ 
sponse (Group et al., 2012; He, Li, Heck, Lustigman, & Jiang, 2006). 
Differences in hydrophobicity and antigenicity were detected in six 
parts in the SI domain (aa 1-552) among the five isolates. The SI 
domain plays an important role in host receptor binding and contains 
a neutralizing epitope (Shang et al., 2018; Woo et al., 2010). These 
results might support the proposition that the chimeric Vietnam 
PDCoV isolates are evolving towards improved host infection and 
immune response evasion. 

To investigate the infectivity of the chimeric Vietnam PDCoV 
(P29_15_VN_1215) isolate compared to the Thai isolate (NT1_1215), 
both virus isolates were propagated in ST cells. The viral growth 


curves were developed and compared to each other. The results 
showed that the growth curve of chimeric Vietnam PDCoV isolate 
(P29_15_VN_1215) in ST cells was relatively lower compared to the 
Thai isolate at 12-60 hpi. This finding suggested that the chimeric 
Vietnam PDCoV isolate (P29_15_VN_1215) had less infectivity in ST 
cells than the Thai isolate. 

To demonstrated whether or not this chimeric Vietnam PDCoV 
(P29_15_VN_1215) isolate can induce clinical disease, the virus 
was orally inoculated in piglets. Clinical diseases were compared 
to piglets that were orally inoculated with the Thai isolate. The re¬ 
sults demonstrated that both PDCoV isolates can induce clinical 
diseases and pathogenicity in piglets. V/C ratio of piglets infected 
with chimeric Vietnam PDCoV was reduced but not significant dif¬ 
ference compared to those of piglets infected with Thai PDCoV. 
The PDCoV infectivity results are in line with other studies. In 
2015, two PDCoV isolates, OH-FD22 and OH-FD100, were also 
shown to cause disease. The infected gnotobiotic pigs showed 
severe diarrhoea, vomiting and severe atrophic enteritis (Jung et 
al., 2015). Moreover, in 2018, piglets were challenged with PDCoV 
KNU16-07 isolate and showed similar clinical signs, namely diar¬ 
rhoea, vomiting and viral enteritis (Jang et al., 2018). The findings 
in our study confirmed that the chimeric Vietnam PDCoV infection 
could induce small intestinal enteritis. 

In conclusion, Vietnam PDCoV is genetically diverse influencing 
from the external introduction. Vietnam PDCoV is further evolved 
into 2 separated subgroups including SEA-1 and SEA-2. Vietnam 
PDCoV in SEA-2 was closely related to Thai and Lao PDCoV. 
Recombination analysis demonstrated that isolates in SEA-1 were a 
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chimeric virus. P12_14_VN_0814, the first Vietnam isolate, and US 
PDCoV isolates were major and minor parents, respectively. The chi¬ 
meric Vietnam PDCoV could induce clinical diseases. 
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